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ABSTRACT 


In  airborne  An t i -Su bma r i ne  Warfare  operations  there  is 
i 
a  critical  requirement  for  maintaining  an  accurate  relative 

plot  of  the  sonobuoys  with  respect  to  the  aircraft.  This 
study  proposed  a  method  for  locating  sonobuoys  in  a  pat- 
tern using  aircraft-to-buoy  slant  range  information.  The 
method  did  not  use  t r i angu 1  a t i on  procedures  and  attempted 
to  minimize  the  restrictions  placed  on  the  aircraft.  The 
study  showed  the  feasibility  of  the  proposed  methodology  and 
the  approximate  errors  to  be  encountered. 
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I .   INTRODUCTION 


The  fixed  wing  aircraft  plays  an  important  role  in  air- 
borne antisubmarine  warfare  (ASW)  today.   The  primary  tacti- 
cal employment  of  the  aircraft  is  through  the  use  of  active 
and  passive  electronic  sensors,  the  most  useful  of  which  is 
the  sonobuoy.   Sonobuoys  are  placed  by  the  aircraft  in  vari- 
ous geometric  patterns  in  the  water  as  the  tactical  situa- 
tion predicts.   The  primary  phases  in  the  prosecution  of  a 
contact  are  search,  detection,  classification,  localization 
and  attack.   Except  in  the  rare  case  of  a  visual  sighting  of 
a  submarine,  the  first  four  phases  are    almost  exclusively 
conducted  by  the  use  of  passive  and  active  sonobuoys  which 
transmit  data,  analyzed  by  sensor  equipment  on  board  the 
aircraft.   This  information  is  then  synthesized  by  the  air- 
craft tactical  commander  to  yield  an  assessment  of  the  tac- 
tical situation. 

One  major  problem  in  airborne  ASW  is  the  lack  of  accu- 
rate a i r craf t/sonobuoy  referencing.   Surface  wind,  sea  cur- 
rents and  errors  in  navigation  quickly  disperse  the  actual 
buoy  pattern  location  from  its  assumed  location.   Moreover 
the  buoys  within  the  pattern  frequently  drift  in  a  random 
manner,  particularly  when  the  distance  between  buoys  is 
great.   It  is  necessary  to  periodically  update  the  actual 
buoy  positions  by  a  visual  or  electronic  "mark  on  top." 
The  flexibility  of  the  aircraft  is  obviously  reduced  by  this 
requirement.   In  the  fine  localization  phase  where  a  combi- 
nation of  two  or  more  buoys  in  close  proximity  to  the  target 


is  providing  the  aircraft  with  fix  data,  it  is  necessary 
to  fly  outbound  to  the  fix   after  passing  over  the  closest 
buoy  providing  data.   This  tactic  is  necessary  to  provide 
the  aircraft  with  final,  accurate  buoy  location  information 
just  prior  to  the  attack.   Frequent  missed  kills  are  attri- 
buted to  attacking  a  fix  position  generated  by  buoys  whose 
accurate  locations  with  respect  to  the  aircraft  were  not 
known  or  recently  updated. 

A  better  method  for  determining  the  location  of  sono- 
huoys  is  required.   A  method  for  locating  buoy  coordinates 
relative  to  the  aircraft,  remote  from  the  buoy,  would  elimi- 
nate the  need  for  "on  top"  repositioning.   Such  a  method 
would  also  offer  the  additional  aircraft  flexibility  re- 
quired for  the  eventual  increase  in  the  tactical  load  placed 
on  the  crew. 

A  system  under  development  [Ref.  1,2]  will  provide  for 
a  slant  range  measurement  to  the  buoy.   This  slant  range 
would  be  more  accurate  than  the  aircraft  navigation  system 
and  could  be  converted  to  horizontal  range  knowing  the  air- 
craft altitude.   The  range  to  a  buoy  can  be  taken  at  any 
time  after  it  has  been  placed  in  the  water  up  to  the  life 
of  the  buoy  and  in  any  tactical  environment  providing  the 
aircraft  is  at  sufficient  altitude  to  receive  the  signal. 
The  use  of  the  ranges  provided  by  the  buoys  under  the  pro- 
posed system   [Ref.  1]  serves  as  the  subject  of  this  study. 

While  t r i angu 1  a t i on  methods  obviously  provide  a  solu- 
tion to  the  buoy  location  problem  [Fig.  1]  th   aircraft 
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A/C  Track 


R.   =   horizontal  range  from  A/C  to  buoy 


A 
B 


A/C  location  on  track  when  first  range 
measu  red 

A/C  location  on  track  when  second  range 
mea  su  red 


0    =   separation  angle  between  range  measurements 


FIGURE    1 
STANDARD  TR I ANGU LAT I  ON  PROCEDURE 
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track,  line  AB,  between  the  buoy  rangings  required  to  com- 
plete t r i angu 1  a t i on  should  be  approximately  equal  to  the 
range  of  the  buoy.   The  buoy  fix  position  is  determined  by 
the  intersection  of  the  range  arcs  to  the  buoy  from  the 
aircraft  position.   The  errors  realized  in  the  buoy  fix 
position  are,  for  the  most  part,  due  to  inaccuracies  in  the 
ranging  system  and  inaccuracies  in  navigation  from  points 
A  to  B.   If  the  leg  AB  [Fig.  1]  representing  the  aircraft 
track  is  allowed  to  be  in  the  order  of  length  equal  to  the 
range  to  the  buoy  such  that  the  range  arcs  intersect  as 
close  to  perpendicular  as  possible,  then  the  errors  due  to 
navigation  and  range  measurement  would  be  transferred  to  the 
buoy  fix  error  in  a  nearly  one-to-one  ratio.   If,  however, 
the  aircraft  track,  leg  AB ,  between  range  measurements  is 
short,  inferring  an  angle  6  of  only  a  few  degrees,  then  the 
range  arcs  will  intersect  nearly  parallel.   The  bearing 
error  of  the  resulting  buoy  fix  position  may  be  extremely 
1 arge . 

tn  that  the  length  of  the  aircraft  track  required  for 
optimal  error  reduction  in  t r i angu 1  a t i on  necessarily  re- 
stricts the  flexibility  of  the  aircraft,  this  study  will 
focus  on  a  method  of  minimizing  the  error  by  statistical 
methods  using  a  rapid  series  of  rangings  to  the  buoy  over 
a  short  interval.   It  is  felt  that  the  statistical  methods 
of  error  reduction  will  allow  an  improved  solution  to  the 
problem  over  standard  tri angulation. 
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f I .   NATURE  OF  THE  PROBLEM 

The  major  portion  of  this  study  is  centered  around 
(1)  determining  a  method  for  locating  sonobuoys  relative  to 
the  aircraft,  (2)  feasibility  of  the  system  with  respect  to 
the  deployment  of  the  aircraft  in  a  tactical  environment, 
and  (3)  the  errors  anticipated  and  how  sensitive  the  pro- 
cedure is  to  these  errors.   A  residual  effort  conducted  in 
conjunction  with  the  development  of  the  sonobuoy  location 
problem   provides  a  conceptual  solution  method  for  deter- 
mining an  accurate  aircraft  ground  speed  and  true  course. 
The  development  concepts  for  the  ground  speed  and  true 
course  determination  are  not  complicated  and  should  therefore 
be  adapted  easily  to  any  aircraft  in  the  fleet  today  which 
would  carry  the  associated  sonobuoy  ranging  equipment. 

Range  information  from  the  buoys  will  be  available  con- 
tinuously, sampled  as  desired.   The  slant  range  will  be 
converted  to  horizontal  range  as  described  by  Ellis  [Ref.  1, 
Fig.  k ] .   There  is  no  inherent  ambiguity  problem  such  as 
mirror  images  of  fix  positions  and  the  solution  is  not 
quadrant  or  half  plane  restricted.   Minor  restrictions  are 
placed  on  the  aircraft  for  short  periods  of  time.   The  air- 
craft, in  most  cases,  will  be  required  to  fly  straight  and 
level  at  a  con.-.tant  true  airspeed  while  sampling  ranges  to 
the  buoys;  however  these  intervals  are  qu'rte  short,  four  to 
twenty  seconds  as  the  tactical  situation  dictates.   This 
restriction  is  not  considered  to  be  excessive  as  will  be 
explained  in  the  formal  development  of  techniques  for 
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sonobuoy  location  and  ground  speed  determination.   In  cer*- 
tain  cases  the  measured  ranges  to  the  sonobuoys  will  be  the 
only  input  required  by  the  system,  which  generates  a  very 
appealing  aspect  to  the  solution  technique.   In  other  cases 
the  aircraft  coordinates  as  well  as  the  measured  range  in- 
formation will  be  required  as  the  only  inputs. 

It  should  be  noted  here  that  navigational  errors  pro- 
ducing erroneous  geographical  track  of  the  aircraft  are  not 
critical  in  the  problem  solution  as  such.   The  primary  pro- 
blem being  attacked  by  this  study  is  that  of  locating  the 
sonobuoy  relative  to  the  aircraft.   If,  in  fact,  the  DR 
track  is  erroneous  then  the  positions  determined  for  the 
buoys  will  be  erroneous;  however,  their  positions  relative 
to  the  aircraft  will  be  accurate  as  they  will  be  plotted 
relative  to  the  aircraft  DR  position.   Accurate  geographical 
track  and  sonobuoy  location  information  can  be  updated  sub- 
sequent to  the  attack  phase  of  the  problem.   Experience  has 
shown  that  two  aircraft  in  the  process  of  relieving  each 
other  may  have  a  relative  error  of  up  to  15  miles  in  geo- 
graphical DR  tracks.   Therefore  the  transfer  of  geographical 
positions  of  buoys  within  a  pattern,  is  frequently  unsatis- 
factory.  When  it  is  necessary  to  turn  the  problem  over  to 
a  relieving  aircraft,  it  will  be  easier  to  transfer  a  rela- 
tive picture,  that  is,  range  and  bearing  of  all  other  buoys 
in  the  pattern  from  one  reference  buoy  close  to  the  relieving 
aircraft. 

Although  the  P  3  -  C  aircraft  will  far  surpass- any  existing 
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ASW  airborne  weapon  system  in  the  fleet  today  it  is  not 
anticipated  that  the  models  used  currently  will  be  phased 
out  in  the  near  future.   It  is  important  to  consider  devel- 
oping a  system  which  may,  with  minimal  modification,  be 
adapted  to  current  models  of  ASW  aircraft.   The  ground  speed 
and  true  course  determination  is  developed  with  this  in  mind 
Considering  the  contribution  of  locating  sonobuoys  from  re- 
mote positions  to  the  prosecution  of  a  contact,  it  may  be 
well  worth  while  to  consider  the  installation  of  a  limited 
purpose  computer  in  the  current  model  VP  aircraft  vice  more 
costly  alternatives. 
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III.   BUOY  RANGING  CONCEPT 

The  sonobuoy  has  a  transmitter  built  in  for  the  purpose 
of  transmitting  to  the  aircraft  the  sensor  signals  received 
from  the  sonobuoy  hydrophones.   The  carrier  frequency  asso- 
ciated with  this  signal  is  channelized  in  each  buoy  for 
buoy  discrimination.   A  buoy  ranging  signal  described  in 
detail  by  Ellis  [Ref.  1]  will  be  transmitted  over  the  exist- 
ing aircraft  receivers.   No  new  receivers  or  transmitters 
for  the  signal  will  be  required.   A  modification  to  the 
detector  section  of  the  sonobuoy  receiver  will  be  required 
for  signal  separation. 

The  sonobuoy  ranging  signal  is  generated  by  a  crystal 
oscillator  placed  in  the  buoy.   The  method  of  range  deter- 
mination is  by  comparison  of  the  phase  of  the  buoy  signal 
with  that  of  a  standard  frequency  generated  on  board  the 
aircraft.   The  phase  difference  is  converted  to  slant  range. 

The  sonobuoy  undergoes  a  warmup  period  in  the  aircraft 
prior  to  drop.   It  generates  a  stable  frequency,  f  ,  and  is 
compard  with  the  aircraft  frequency,  f  ** .   The  phase  dif- 
ference at  this  time  and  at  zero  distance  is  measured,  AG  • 

o 

At  any  time,  t,  after  the  sonobuoy  is  dropped  in  the  "water 

a  new  comparison  is  made  between  f   and  f " ,  and  the  phase 

r       r  r 

AG  ,  is  determined.   The  resulting  phase  difference,  AG  , 

where   AG   =  AG,  -  AG  ,  i 5  then  converted  to  slant  range: 
r      t      o  _ 

subsequently  to  horizontal  range.  The  buoy  range  frequency 
has  a  wavelength  of  approximately  75  kilometers.  One  cycle 
difference  in  phase  of  the  buoy  signal  from  the  "reference 
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is  equivalent  to  75  km.  distance.   A  buoy  at  a  range  of 
10  km.  and  85  km.  would  yield  the  same  phase  difference 
causing  an  ambiguity.   It  must  be  known  what  multiple  of 
75  km.  the  buoy  is  from  the  aircraft.   It  is  not  anticipated 
that  this  ambiguity  will  cause  erroneous  ranging. 

The  range  errors  generated  by  the  system  are  covered 
in  detail  by  Ellis  [Ref.  1]  but  one  error  of  primary  con- 
cern in  the  tactical  employment  of  the  system  is  the  linear 
increase  in  frequency  with  respect  to  time  within  the  buoy 
oscillator.   Complete  frequency  stabilization  is  not  yet  a 
"state-of-the-art"  realization.   Testing  to  date  with  the 
oscillator  currently  in  use  in  the  buoy  has  generated  an 
average  drift  rate  comparable  to  a  range  error  of  10^0  meters 
per  hour.   If  a  means  was  devised  for  the  subtraction  of  the 
linear  drift  rate  as  a  function  of  time,  the  remaining 
errors  in  the  system  could  be  aggregated  to  a  classification 

2 

of  normal  with  a  mean  of  zero  and  a  variance  of  a  .   This 

£ 
2 

aggregation  of  errors  being  N(0,a  )  would  allow  a  broad 
application  of  statistical  methods  for  the  solution  of  the 
error  reduction  problem  in  the  tactical  employment  of  the 
buoy  range.  .  It  is  therefore  assumed  that  in  the  following 
study  the  average  linear  error  rate  is  eliminated  from  the 
problem  and  the  assumption  that  the  remaining  errors  are 

2 

N(0,a  )  will  motivate  the  development. 
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IV  .   SUMMARY  OF  PROBLEM  SOLUTION  TECHNIQUES 

The  following  methods  are  developed  for  the  solution 
of  (l)  Sonobuoy  Location,  (2)  Aircraft  Ground  Speed  and  (3) 
Aircraft  True  Course.   This  section  provides  a  non-techni- 
cal summary  of  the  problem  geometry,  basic  solution  approach 
and  the  assumptions  and  restrictions  governing  the  technique 
used  in  each  of  the  three  problem  areas  mentioned  above.   It 
is  designed  as  a  background  preface  to  the  more  detailed 
mathematical  development  and  should  by  itself  acquaint  the 
reader  with  a  complete  broad  picture  of  what  has  been  accom- 
plished.  The  mathematical  development  will  be  included  in 
a  subsequent  section. 

A.    SONOBUOY  LOCATION 

It  is  desired  to  provide  a  system  whereby  the  aircraft 
may,  in  any  phase  of  the  -tactical  prosecution  of  a  contact, 
locate  any  or  all  of  the  sonobuoys  in  the  pattern  with  re- 
spect to  the  aircraft.   To  provide  the  aircraft  with  the 
flexibility  needed  in  the  more  acute  phases  of  the  problem, 
the  system  must  not  place  heavy  restrictions  on  the  maneuver 
ability  of  the  aircraft  such  as  is  required  in  standard 
t r i angu 1  a t i on  techniques.   A  method  whereby  the  aircraft 
could  obtain  locations  from  a  single  position  and  time  or 
with  a  minimum  advance  leg  would  certainly  be  preferred. 

Recalling  the  discussion  in  the  1  a s t  ~pa rag ra ph  of  the 
introduction,  an  extended  aircraft  track,  line  AB  [Fig.  1], 
would  have  to  be  flown  in  order  that  9  be  sufficiently  large 
to  reduce  fix  error  of  the  buoy  position.   In  a  tactical 
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situation  it  may  not  be  feasible  for  the  aircraft  to  fly 
the  leg  AB .   If  the  aircraft  desired  to  maneuver  between 
poiiits  A  and  B  it  must  do  so  at  the  possible  cost  of  addi- 
tional navigational  errors  as  aircraft  navigation  becomes 
degraded  in  many  cases  as  a  result  of  maneuvering. 

The  geometry  in  Figure  2  was  considered  as  a  starting 
point  for  the  development  for  sonobuoy  location.   This  basic 
geometry  will  prevail  in  the  description  of  the  problem. 
Consider  an  aircraft  with  a  buoy  in  the  relative  position 
shown.   The  aircraft  is  required  to  hold  any  true  heading 
and  airspeed  for  as  long  as  it  is  required  to  measure 
ranges  Ri,  R2  and  R3  [Fig.  3].   The  time  interval  between 
each  range  measurement  is  short,  approximately  five  seconds. 
At  200  kts.  an  aircraft  travels  111  yds. /second.   As  many 
ranges  as  desired  may  be  measured  but  more  time  is  required 
and  the  additional  accuracy  obtained  by  more  ranging  is  not 
significant  as  a  trade  off  in  a  critical  tactical  situation. 
This  will  be  a  tactical  decision.   The  study  uses  three 
range  measurements  throughout  the  descriptive  phase.   Addi- 
tional ranges  were  used  in  testing.   It  is  assumed  that  in 
the  short  period  of  time  the  aircraft  is  required  to  Tly 
straight  and  at  a  constant  air  speed  that  the  coordinates 
(A.,B.)  [Fig.  31  are  accurate  relative  to  each  other.   Any 
inaccuracies  occuring  will  be  absorbed  in  final  fix  error. 

The  ranges  R.  and  the  aircraft  coordinates  (A.,B.)  are  in- 
1  11 

puts  and  are  assumed  to  be  known.   The  buoy  coordinates 
(X,Y)  are  unknown. 
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CPA 


-4- 


A/C  Heading 


R.   =   horizontal  range  from  A/C  to  buoy 

CPA  =   closest  point  of  approach  to  buoy 
with  respect  to  the  A/C 


relative  bearing  of  buoy  with  respect 
to  the  A/C 


FIGURE   2 
BASIC  TWO  DIMENSIONAL  PROBLEM  GEOMETRY 
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(A1  .B,) 


(A2,B2) 


(A3,B3) 


(X,Y) 


Head  i  ng 


R.       =  horizontal  range  to  buoy  as  measured  by  A/C 

1  =  1  ,2,.  .  .  ,n 

(A .  ,  B . )  =  A/C  coordinates  at  time  of  measurement  of 
i   i 

range  R.  to  buoy  .  ,  9 

(X,Y)    =  coordinated  of  buoy  (unknown) 


FIGURE   3 
S0N0BU0Y  LOCATION  PROBLEM  GEOMETRY 
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The  model  used  to  determine  the  (X,Y)  position  of  the 

buoy  is  the  method  of  Least  Squares.   The  basic  assumption 

I  2 

forithe  use  of  this  model  is  that  the  errors  are  N  ( 0  ,  a  ). 

The  approach  is  to  minimize  the  sum  of  the  squares  of  the 

errors  between  the  measured  ranges  to  the  buoy  and  the 

actual  ranges.   The  equation  for  the  Sum  of  Squares; 


SS   =   J  (O.-E.) 
1-1   '    ' 


i  =  1 


where:     O.=observed  value 
i 

E.=expected  value 
i 

0.-E.=e.  (error  term) 
i   i   i 


is  adapted  to  the  sonobuoy  ranging  problem; 


SS  =  I    (r.-rD 

i  =  1   '   ' 


where;     R.=measured  range 

R.=expected  value 

=  [(X-A.)2H-(Y-B.)2]i 

R. -rT=e. 
i   i   i 


The  approach  of  minimizing  the  sum  of  the  error  squared 
termsis  that  of  evaluating  the  first  derivative  of  SS  at 
zero  with  respect  to  the  X  and  Y  values  desired; 


SS 


I  (e.) 


4  2 


i  =  1 


I    {R,  -  [(X-A.)   +  (Y-B.)  ]*} 
1=1   '  '  ' 


ass 


8SS 
8Y 


=  0 


This  development  produces  a  set  of  normal  equations  from 
which  the  values  of  X  and  Y  are  determined  by  an  iterative 
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p  roces  s  . 

The  iterative  process  requires  that  an  estimate  of  the 
buoy  angle  relative  to  the  aircraft,  6,  be  supplied.   This 
initial  angle,  9o,  provides  a  starting  point  for  the  pro- 
cess which  by  making  alternately  successive  approximations 
of  X,  Y  and  0  the  process  converges  on  the  values  of  X  and 
Y  with  the  minimum  error  possible  for  the  geometry.   Values 
of  60  up  to  and  including  ±  15   in  error  of  the  actual 
angle  9  were  tried  [Fig.  k]  .       There  was  little  effect  on 
the  rapid  convergence  of  X  and  Y  to  within  acceptable  values. 
In  consideration  of  the  latitude  in  choosing  a  value  of  0o 
it  is  felt  that  the  bearing  supplied  to  the  buoy  by  the  DME 
indicator  is  satisfactory  as  an  input  value.   The  value  of 
So  was  considered  to  have  been  chosen  as  an  estimate  of  the 
0  for  the  mi:' -range  measurement. 

Two  studies  were  conducted  to  determine  the  feasibility 
of  the  Least  Squares  model  for  the  solution  of  the  sonobuoy 
location  problem.   The  first  study  was  a  sensitivity  analysis 
to  determine  the  effect  of  the  following  on  the  convergence 
rate  of  the  process;  (1)  distance  to  the  buoy,  (2)  relative 
angle  of  the  buoy,  0,  (3)  values  of  6d  as  an  estimate  of 
and  (h)     number  of  ranges  used  in  the  solution  of  the  pro- 
blem. 

The  second  study  was  an  error  analysis  to  examine  the 
effect  of  range  error  on  the  final  values  of  buoy  coordi- 
nates, X  and  Y.   Errors  of  varying  magnitudes  were  applied 
to  the  ranges  in  different  combinations  at  different  relative 
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Head  i  ng 


A   =   A/C  location  when  horizontal  range  to  buoy 
measu  red  .  - 


9   =   actual  relative  angle  of  buoy  with  respect 
to  A/C  heading 


8   = 
o 


estimated  values  of  relative  angle  of  buoy; 
error  applied  in  1°  increments  up  to  +  15° 
of  actual  relative  angle 


FIGURE  k 

RELATIONSHIP  BETWEEN  ESTIMATED  AND  ACTUAL  VALUES  OF  RELATIVE 
BUOY  ANGLE  WITH  RESPECT  TO  A/C  HEADING 
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buoy  angles.   The  results  of  these  studies  are  explained  in 
detail  in  Appendix  A. 

B.    GROUND  SPEED  DETERMINATION 

The  development  of  the  method  of  ground  speed  deter- 
mination is  based  on  the  concept  that  if  an  aircraft  was 
approaching  a  buoy  (0  relative  bearing)  which  was  providing 
accurate  range  information  then  the  change  in  the  horizontal 
range  to  the  buoy  with  respect  to  time  would  be,  exactly, 
the  ground  speed.   This  is  a  fast,  accurate  and  simple  way 
to  determine  ground  speed  if  it  could  be  ascertained  that 
the  relative  bearing  was  exactly  6   relative.   This  is 
difficult  to  do  as  the  DME  indicator  does  not  provide  that 
accuracy.   As  the  relative  angle  of  the  buoy  increases  from 
zero  the  first  derivative  of  the  range  to  the  buoy  with  time 
becomes  a  function  of  the"  cosine  of  the  relative  angle.   If 
there  existed  a  scheme  for  determining  the  change  in  range 
to  the  buoy  with  respect  to  time  with  the  relationship  of 
the  relative  angle  then  the  ground  speed  could  be  determined 
for  any  relative  angle  for  the  buoy. 

In  Figure  5  a  geometry  similar  to  that  of  the  previous 
section  is  used.   Consider  an  aircraft  approaching  a  buoy 
as  indicated.   Were  the  aircraft  to  continue  on  this  head- 
ing it  would  pass  the  buoy  at  its  closest  point  of  approach 
(CPA)  at  a  point  perpendicular  to  the  air c> aft  track.   In 
the  problem  solution  the  aircraft  chooses  any  buoy  whose 
relative  bearing  is  close  to  the  aircraft  heading.   The  air- 
craft is  then  required  to  fly  straight  and  level  at  a 
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Y  =  CPA 


CPA   =   closest  point  of  approach  to  buoy  with 
respect  to  the  A/C  heading 

R.    =   horizontal  ranges  to  the  buoy  measured 
at  A/C  locations  1,  2  and  3 

C     =   equal  distances  (assumed)  generated  by 
flying  equal  times  at  a  constant  true 
airspeed  between  points  1,  2  and  3 


Note:    Buoy  coordinates  are  not  known,  moreover 
are  not  required  for  this  problem. 


FIGURE   5 
GROUND  SPEED  DETERMINATION  PROBLEM  GEOMETRY 
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constant  true  air  speed  for  a  time  interval  sufficient  to 
measure  from  three  to  eight  ranges  to  the  buoy  in  question. 
The  time  between  each  measurement  is  short,  two  to  four 
seconds.   The  time  between  each  range,  C,  should  be  exactly 
equal.   This  could  be  done  by  electronic  time  gating.   This 
time,  C,  also  represents  equal  distance  between  range  mea- 
surements.  It  is  assumed  that  for  short  periods  of  time  and 
by  electronic  time  gating,  these  distances  will  be  equal. 
Any  errors  are  absorbed  in  the  total  ground  speed  error. 
The  actual  or  estimated  position  of  the  buoy  is  not 
required  as  an  input  to  the  problem.   Also  it  is  not  re- 
quired to  know  the  heading  or  airspeed  of  the  aircraft, 
only  that  they  are  constant.   The  evaluation  of  C  deter- 
mines the  ground  speed.   The  use  of  the  Pythagorean  theorem 
produces  three  equations  for  the  three  ranges  measured; 


/      2     2     2 

(X+2C)   +  Y   =  Ri     for  range  measured  at  1 


(X+  C)2  +  Y2  =  R2 


for  range  measured  at  2 


2         22 
X       +  Y   =  R3     for  range  measured  at  3 


The  magnitude  of  the  ranges  only  are  known  and  not  the  d i rec 
tion.   There  are  three  equations  and  three  unknowns;  X,  Y 
and  C.   The  values  for  C  are  evaluated  from  the  above  equa- 
tions for  the  three  ranges; 


C   = 


2     22 
Ri-  2R2+  R  3 


By  a  similar  development  C  is  determined  for  any-number  of 
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ranges  measured.   Up  to  eight  ranges  are  used  in  the  de- 
tailed mathematical  development  included  later. 

I  Two  studies  were  conducted  to  determine  the  feasibility 
of  the  ground  speed  determination  method.   The  first  study 
was  a  sensitivity  analysis  to  determine  the  effect  on  ground 
speed  error  by  varying  the  following;  (1)  number  of  ranges 
used,  (2)  time  between  range  measurement,  (3)  ranges  to  the 
buoy  and  (4)  the  relative  bearing  of  the  buoy. 

The  second  study  was  an  error  analysis  on  the  ground 
speed  as  affected  by  errors  in  the  measured  ranges.   The  re- 
sults of  these  studies  are    combined  in  Appendix  B. 

This  method  for  determining  ground  speed  is  expected  to 
have  limited  application,  however,  because  of  the  simplicity 
of  concept  it  should  be  quite  reliable  to  supplement  break- 
down of  more  sophisticated  inei  tial  and  dopplar  navigational 
systems.   Frequently,  in  fine  localization,  a  timed  run  out 
from  some  reference  to  a  fix  position  for  attack  is  required. 
It  is  felt  that  this  ground  speed  determination  method  can 
be  utilized  on  the  final  leg  to  more  accurately  determine 
the  ground  speed  and  hence  a  more  accurate  time  to  fix  on 
that  leg.   The  aircraft  will  maintain  a  true  course  and  air- 
speed to  that  fix  position  which  meets  the  requirements  for 
the  solution. 
C.    TRUE  C0URSE/  DETERMINATION 

It  was  felt  that  the  use  of  the  range  measurement  to  a 
buoy  as  the  only  input  could  provide  additional  information, 
if  applied  properly,  in  evaluating  a  "drift  vector".   Consider 
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in  Figure  6-a,  an  aircraft  flying  a  true  heading  of  due 

North.   This  true  heading  is  held,  at  a  constant  true  air 

speed  (TAS) ,  for  a  time  t w .   Vector  D.,  represents  the  true 

air  speed  times  the  time  flown,  t...   The  aircraft  then 

turns  to  a  true  heading  of  East  and  holds  this  heading,  at 

the  same  TAS,  for  a  time  t  _.  where  tr  =  t .,  .   While  on  each 

E         E     N 

of  these  headings  the  aircraft  determines  ground  speed  by 
methods  described  in  the  previous  section.   A  different 
buoy  is  used  for  ground  speed  determination  for  each  of  the 
two  headings.   If  a  wind  from  the  Southeast  were  affecting 
flight,  the  ground  speed  for  the  northerly  leg  would  be 
greater  than  the  TAS  and  the  ground  speed  for  the  easterly 
leg  would  be  less  than  the  TAS,  [Fig.  6-b].   The  triangle 
Qi    VI    D^.    shows  the  second  leg  geometry  without  the  advance 
and  transfer  of  the  aircraft  in  the  turn.   The  magnitude 
only  and  not  the  direction  of  these  ground  speed  vectors  is 
known.   It  should  be  noted  that  any  initial  true  heading 
and  TAS  could  be  flown  and  a  subsequent  right  or  left  turn 
could  be  executed.   Maximum  accuracy  of  the  drift  vector  is 
obtained  when  the  subsequent  heading  is  perpendicular  to  the 
initial  heading. 

In  the  development  for  the  determination  of  a  drift 
vector,  the  advance  and  transfer  of  the  aircraft  was  elimi- 
nated in  further  diagrams  and  does  not  alter  the  geometric 
accuracy  of  the  development.   The  actual  path  of  the  air- 
craft, less  the  advance  and  transfer,  [Fig.  7~a],  is  repre- 
sented by  the  vectors  G..  and  Gr  whereas  the  vectors  D..  and 

N       t  N 
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D. 


Advance  and  transfer 
of  the  A/C  in  turn 


Fig.  6-a 


Wind 


Fig.  6-b 


D..  =  t 


Dr  =  t 


VL  =  t, 


W„  =  t 


G.,  =  t 


GP  =  t 


TAS 


TAS 


where  t..  =  t  _. 
N     E 


DN  "  DE 


W,  =  w2 


(Windspeed  +  additional  errors) 

(Windspeed  +  additional  errors) 

(Groundspeed  for  North  leg)  magnitude  only 

(Groundspeed  for  East  leg)  magnitude  only 


FIGURE   6 
TRUE  COURSE  DETERMINATION  PROBLEM  GEOMETRY 
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(X,Y) 


Fig.  7~a 


(c,d) 


(e.f) 


(a,b) 


W  =  drift  vector  =  wind  vector  plus  additional  error 

components 

Known  Values : 

A/C  assumed  position  coordinates,  (a,b),  (c,d)  and  (e,f) 
vectors  D^  and  D^  ,  both  magnitude  and  direction 
vectors  G^  and  G ..  ,  magnitude  only 

Unknown  Values: 

vectors  G^  and  Gr,  direction  only 

coord  inates  (X,Y) 

drift  vector,  W,  both  magnitude  and  d  t-rect  ion 


FIGURE   7 
TRUE  COURSE  DETERMINATION  SOLUTION  GEOMETRY 
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D   represent  the  "no  wind"  track.   The  geometry  of  the  pro- 
blem was  altered  for  the  purpose  of  the  mathematical  solu- 
tion technique.   The  triangle  GpW„D   was  shifted  [Fig.  7~b] 
such  that  W  _  coincided  with  W.  since  W  _  =  W  -  .   The  vector 
D_  did  not  change  in  orientation  thus  the  relationship  be- 
tween the  two  triangles,  GrW_Dr  and  G..W.D..,  was  not  changed 

E  2  E       N  1  N 

The  resulting  form  revealed  that  the  vectors  G..  and  Gr  had 
3  N       E 

the  same  terminal  point.   The  known  and  unknown  values  in 
the  figure  are  listed  at  the  bottom  of  Figure  7- 

A  pair  of  quadratic  equations  were  used  to  solve  for 
the  terminal  point,  (X,Y); 


(X-a)2+  (Y-b)2  =  G2 
(X-c)2+  (Y-d)2  =  G2 


The  magnitude  of  the  d  r  i  f~t  vector,  W,  was  then  solved,  as 
with  the  determination  of  the  coordinates  (X,  Y)  both  the 
coordinates  of  its  origin  and  terminal  point  were  known. 


W   =  ^  (X-e)2+  (Y-f) 


The  resulting  drift  vector  was  then  used  in  conjunction  with 
TAS  and  true  heading  to  determine  an  actual  course  and 
g  round  s  peed  . 
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V.   SONOBUOY  LOCATION  TECHNIQUES 

From  the  description  of  the  basic  sum  of  squares  model 
in  section  IV,  the  sum  of  squares  equation  for  the  parti- 
cular geometry  in  Figure  3  was  obtained; 


1 


2   i   2 


SS  =  ±      I    {R.-[(X-A.)  +(Y-B  )  ]*} 
i  =  1 


In  further  development  let 


x.  =  X-A. 
i        i 


y.  =  Y-B 


(1) 
(2) 


then 


1 


2    2   i   2 


SS  =  -  I     [R.-(x.+y.)2] 
i  =1 


(3) 


Equation  (3)  was  then  differentiated  with  respect  to  X  and 

v,    ,     i4.j  •      9SS  3  S  S    _ 

Y  and  evaluated  at  zero,  i.e.,   ~ v  =  ~ 5  =  0 

o  A  d  Y 


in  =  l.VRi-^i+v;)i](i)(2x.)(x%y;)-i  =  o 

i  =  1 


n  i-i 


R.x 


2    2 
(x.+y.)^ 


1     °   j=1   ' 


=  0 


W 


Let  (a)  b< 


,   n     R.x. 
—   Y 

n    L 


2    2   i 

'"1  (x.+y.)* 
I  '   i 
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and  let  (b)  be 


I      1 


n  . L 


x  . 


i  =  1 


Each  term  of  the  expression  (a)  was  multiplied,  both  numera- 


tor and  denominator,  by  — 

'  x  . 


1  I 

n  .  L 


R.x.  (7  ) 
1  1  x  - 


1 


1 


1     221    n  .  , 


u  4      2      n  .  L 


x . 

1 


i  =  1  (1+tan20.)2 


where  6.  is  the  relative  angle  of  the  buoy  at  the  time  the 
range  R.  was  measured,  and 


y.  Y-B. 

tan  6 .  =  —  =  v — T 

1    x .  X-A . 

t  1 


(5) 


Thus  (a )  becomes  ; 


1  1 

n  .  L 


R. 


i=l  (Sec  0.  ) 


1   n 

—  =  -  I      R.co 

-i-  n     «       1 


n  .  " ,   1      1 
1  =  1 


The  second  term,  (b)  was  expanded; 


1   ?      1   n 


X   ?  "   1 


j-    ;x.si    1   (x-a.)=  ^-    y   -  L  I   A.=  X-  -'-  V  A. 

i=1        i=1  i=1    i=1  i=1 


The  two  terms  (a)  and  (b) ,  were  then  combined; 

n  n 

-    y  R.cose.-x  +  -    y   A.  =  0 

n  .  L .     1      1  n  .  L .     1 

:1  i  =  1 

1   n  n 

-    y  r.cosg.  -  -  y  A. 


i=1 


=  1 


i  =  1 


(6) 


3* 


ass 


A  similar  development  for    v  =  0  yielded 


n  1   n 

1  I    R.sinO.  -Mb. 

i=l  i=1 


(7) 


The  unknowns  X  and  Y  could  not  be  solved  directly  as  8 .  is 
a  function  of  X  and  Y  (see  equation  5);  therefore,  an  initial 
value  common  for  all  the  6.  was  estimated,  namely  0  ,  and 
initial  values  for  X  and  Y,  i.e.,  XQ  and  YQ  were  determined 
using  6  for  each  6.  in  equations  6  and  7  •   These  values  of 
X.  and  Y   Ln  turn  were  used  to  update  the  values  for  0., 


name  1 y ; 


tanr.  = 


Yi-»i 

i    X0  -A. 


■V; 


=  1  ,2 ,  .  .  .  ,n 


These  values  for  the  0.  were  then  substituted  in  the  equa- 

i 

tions  6  and  7  to  determine  X   and  Y,.   This  process  was 
repeated  until  the  convergence  of  X  and  Y  to  within  accept- 
able values  was  attained. 

The  foregoing  iterative  process  proved  to  be  lengthy. 
The  angles  0.  were  quite  insensitive  to  small  changes  in 
the  positions  of  X  and  Y  through  the  iterative  process.   It 
was  felt  that  an  application  of  Newton's  method  in  two  vari- 
ables might  provide  a  more  rapid  convergence.   The  method 
involved  the  expansion  of  the  functions  (as  determined  by 
equations  6  and  7  using  0  )  in  Taylor's  Series  in  two 
variables  about  X0  and  YQ.   The  expansion  would  terminate 
after  the  linear  terms  and  the  remainder  ignored.   The 
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The  resulting  equations  would  be  solved  for  X  and  Y  and  the 
solutions  used  for  an  updated  expansion.   This  process  would 
be  repeated  until  the  values  of  X  and  Y  were  within  accept- 
able limits. 

The  sine  and  cosine  of  the  relative  angles  9.  as  ap- 
pearing in  equations  6  and  7  were  developed  in  terms  of  x. 
and  y.  (see  equations  1  and  2); 


cos 


1 


1 


1 


(sec  G.)2 


(1  +  tanV)*    (l+y?/x!)i 

I  '  I     ! 


I 


(x?+y!)2 


sin  0.  = 


1  o  9     9    1 


(esc  G.)2    (1+cot  0.)2   (l+x./y.) 


I 

— 2 — r~ r 

(x.+y.  * 

i   i 


The  equations  6  and  7  were  then  cast  into  the  homogeneous 
forms  be  1 ow ; 


9l(X,Y)  =  X  -  i  I    {R 


%r-  A.}  =  0 


1  =  1   '  (xt+yt) 

i  '  i 


g2(x,Y)  -Y-i.I  ^i-Ai  -  B.}  =  0 


■=1    (x.+y.) 


The  form  of  the  Taylor  expansion  is 


g(X,Y)=g(X0,Yo)+gv(X0  ,Y0)  (X-X Q ) +g  (X q  ,  Y 0 )  (Y-Y 0 ) +Rema i nde r 

x  y 
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Specifically,  g( expanded; 


i         x .  n        y . 

l(X,Y)  =  X0--  I    {R.  r-! r  "A.}+[1--   £  R 


1 


'  r..z...2\i  i 


i=1     Cx.+y . ) 
i  'i 


— ] (X  X  o) 
i  =  1   (.x  .  +y  .  ) 


+  [-  I  R 


X;Y 


'  (.,  2^. .  ^  3/2 


i=1   lx .+y  .  ) 


](Y-Y0)  =  0      (8) 


and  g  2  expanded  ; 


x.y . 


92(X'Y)=Yo-n-  I    {Ri,  2   Tri    "V^n-  .1  Rt,  2  2  8yJ(X-Xo) 


o  -n  .  «■-.,    I  i     2   2 

=1    Cx.+y. J 


2   2.3/, 

i  =  1   (x.+y.)  /2 


+  H-I  j/i-ra-KY-Y.)  -  o       (9) 

1=1   (x.+y.) 


with  x.  and  y.  evaluated  at  X0,  Y0. 

Equating  (8)  and  (9)  to  zero  yielded  the  following  system 

of  linear  equations; 

2 

1   n       Y  •  1   n       x .  y .    v 


1=1   Cx.+y . ) 


n  •  1  I '/   a.  2\  V2 
i=1   (x.+y . ) 


.  y. 


Y0  D 


X0  n 


_  _u  y  R  LL 

n     A     J     2     2   3 A 

(x.+y.)  /2 


n  J/i  J  *   2  £  "j 

i=1   Cx.+y, )        i=1 


i   i 


«   n    x . 

+  L  I 

n    ^ 


1    n 


i =1 Cx .+y . ]       i=1 
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1 


n 


2 
X  . 


,       n  x  .  y . 

[1      I    R.  '     '        JX    +     [1--      I    R. 

L  n       .  L  .       I    ,      2  2  .    3/2 


i     7  i 


i  =  1        (x.+y  .  ) 


]Y 


i  =  1     '  (xT+yt) 


X.n  x.y.  Ynn 


i  /      2        2V  3/2       n      .',     2.      2*  3/2 

i  =  1  (x . +y . )  ' 


•  =  1        (x.+  Y.) 
i       t 


y 


i  i 


n  y 


2— T—  +    -      I    B- 
n     .,/    2.     2\'  n     .  *■  „     i 


i -1  (x. +y . ) * 


Cramer's  rule  was  applied  to  solve  for  X  and  Y.   The  coef- 
ficient matrix,  H,  of  X  and  Y, 


H    = 


Di  D3 


2  i+ 


1=-     I    R. 

n     . L  ,      i 


i  1 


I   R 


X:Y 


.     i  /    2       2  \  3/9        n       '        i  /    2       2.3/, 
=  1        (x.+y.)  /2  i  =  l        (x.+y.  )  /2 


y 


x.y 


i-  y  r  - 

n     .L  ,      \  ,    2  .     2  v  3/2 
i=1        (x.+y.) 


'-1    I   R 

n  ^ 


n  x 


n      .  L  .       I    /      2  2^3/2 

1=1        (x.+y  .  )   f2 
i     7  i 


and  the  vector  of  constants, 


I  R 


x.y. 


x       n 

-    I  R 


i  .  1 


n  x 


i  .  1 


y  ! yL  y  A 

n     sS"'/    2.     2  \  V2       n     .'"ii    2       i\*/z'  t\.L .  ,    2       2.i     n.^, 
i=1       (x.+y. J  i=1        (.x.+y.)  i  =  1(x.+y-.)5       i  =  1 


i        i 


x.     n  x.y.  Y„     n 


I  I 


-    I  R 


U         r< 

—      I    R 


xl  n         y 

+1     1 

n  ^ 


1       n 

n      L' 


n     .  *• ,     i/2J.2N3/2     n     .  «■ ,     i.223/2     n     .  L  22N,     n.L. 

i=1        (x.+y.)  i=1        (x.+y.)  i=1(x.+y.)2       i=1 


i        i 


i        i 


38 


were  applied  in  the  following  manner; 


X  = 


Y  = 


Dc  •  D,-  D,  •  D, 


Di  •  DL  -  D,  '  D. 


1      6      b      z 


D,  •  D  -  D  '  D0 

1      h  3     2 


(10) 


(11) 


Equations  10  and  11  replaced  equations  6  and  7  in  the  itera 
tive  process.   This  method  proved  highly  successful  as  a 
rapid  convergence  process.   The  results  of  the  analysis  of 
this  process  are  represented  by  the  figures  in  Appendix  A. 

With  the  success  of  the  second  method  an  attempt  was 
made  to  interpret  the  distinctive  improvement.  The  coef- 
ficient matrix,  H,  was  broken  down  and  analyzed; 

2 


H    = 


1 


0 


n     .  L . 


1  /*J.*\  3/2 


i=1     '(x.+y.) 


1  I  * 

n     .  L  , 


x.  y  . 


1/      2,     2\V2 


1=1        (x.+y.) 


1 


I  R 


x .  y  . 

i  '  i 


i=1        (x.+y.) 


2       3/2 


1  " 


n     .",     i  /    2       z. 
i  =  1        (x.+y  .  ) 
i        i 


V2 
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The  elements  of  the  second  term  were  then  factored 


"  I 


Ri     ,   .   *    v 

—2 2-T(sin  V 


1=1  Cx .+y . ) 2 


i=1  (x.+y . ) i 
i   i 


j (sin6. cos0. ) 


n       R .  n     R . 

I   — 1 £-7-(sine  cos9  )  -£      '    (cos  0  ) 

=1  (x'+y!)*      '      '   l-l(x?+v?)i 


i 


i  =  1  (xt  +  yf) 


The  matrix  below  extracted  from  the  second  term  of 
the  coefficient  matrix  has  a  determinant  whose  value  is 
equal  to  zero  for  any  values  of  6.; 


-sin  6 


i 


sin6. cos6 

I 


s i no . cos  .    -cos  0 . 
i     i  i 


If,  in  fact,  there  were  no  errors  between  the  measured 

2       2   1 

ranges  and  the  actual  ranges  then  R.  =  (x.  +  y . ) ?  and  the 

R.  iii 

factored  term,  g ,  would  equal  one.   This  implies 

(x.  +  y.)* 
that  if  there  were  no  error  between  the  measured  and  actual 

ranges,  the  second  term  in  the  matrix  H  would  go  to  zero 

l  eav  i  ng ; 


H  = 


1    0 
0    1 


It  appears  that  the  existing  error  between  the  actual  and 
measured  ranges  causes  variations  in  the  value  of  the  second 
term  in  the  coefficient  matrix  H  and  therefore  provides  the 
driving  force  of  the  convergence  process. 


kO 


V  I  .   GROUND  SPEED  DETERMINATION 

From  the  problem  description  in  Section  IV,  the  mathe- 
matical model  below  for  ground  speed  determination  was  de- 
veloped.  Figure  5  provided  the  geometry  for  the  problem. 
The  particular  geometry  in  the  figure  indicates  three  range 
measurements  taken  at  equal  intervals.   Ground  speed,  using 
from  three  to  eight  ranges  measurements  was  developed  in 
this  study.   The  distances  C-calculated  below  are  inputs  to 
the  ground  speed  determination  equations  18  and  19  at  the 
end  of  this-  section  where  i  =  3  ,  ^  ,  •  .  .  ,  8  represents  the  num- 
ber of  measured  ranges  utilized  in  the  calculation. 

For  the  three  range  geometry  the  known  values  are  R  , 
R„  and  R  .   It  is  also  assumed  that  the  time  between   range 

2  3 

measurements  is  the  same  in  each  case.   The  unknown  values 
are  X,  Y  and  C.   The  three  equations  from  which  C  was  deter 
mined  are  applications  of  the  Pythagorean  theorem. 

(X  +  2C)2  +  Y2  =  Rj  ■*  X2  +  kCX    +  kC2    +  y2  =  R2  (a) 

22222  222 

(X    +       C)       +    Y       =    R0  -*    X       +    2CX    +       C       +    Y       =    R2     (b) 


2 
2 

+    Y       +    R,  -*   X 

3 


2  2 

+    Y       =    R3     (c) 


Equations  (a)  plus  (c)  minus  2(b)  ; 


X  +  ACX  +  hC       +  Y   =  R1 

2  2        2 

X  +  Y   =   R3 

2  .            2       2        2 

-2X  -  ACX  -  2C   -2Y   =-2R, 


2C 


Rj  —  2Ro  ^*  R3 


2 

2         2 

-    2R2    +    R3 

2 

(12) 
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The  development  for  four  ranges; 


(X  +  3C)2  +  Y2  =  Rj  —  X2  +  6CX  +  9C2  +  Y2  =  R2  (a) 

(X  +  2C)2  +  Y2  =  R2  -  X2  +  ^CX  +  kZ2  +  Y2  =  R2  (b) 

(X  +   C)2  +  Y  =  R2  -~  X2  +  2CX  +  C2  +  Y2  =  R3  (c) 

X2  +  Y2  =  R2  -  X2  +  Y2  -  R2  (d) 


Equations  (a)  minus  (b)  minus  (c)  plus  (d); 


2  222 

X  +  6CX  +  9C  +  Y   =  R 

?                  2  2       2 

-X  -  i*CX  -  kC  -  Y   =  R^ 

2                 2  2       2 

-X  -  2CX  -   C  -  Y   =  R3 

2  2       2 

X  +  Y   =  R„ 


hC 


2       2      2       2 

=  R,  _  Rp  "Rj  +  Rl 


2       2       2       2 

Rl   ~  R2   "  R3   +  Rk 


(13) 


The  development  for  five  through  eight  ranges  was  omitted 
but  was  similar  to  the  two  cases  above; 


2         2       2       2       2 

2R2  -  3R2  +  R3  -  R,  +  R5 


8 


dM 


n   r.2  nr,2  r,2  r,2        r,  2        _2 

f2R,  -  2R0  -  R0  +  R,  -  Rr  +  R, 

1         2       3       1+56 


12 


(15) 


,2       2       2        2       2       2       2 

R,  -  R   +  R   -  2R,  +  Rc  -  R   +  R„ 

12       3         4       5       6       7 


12 


(16) 


C   = 

e 


T 


[R!  -  R0  +  r!  -  R2  -  R2  +  R2  -  R2  +  R2 

1       2       3       h  5       6       7       8 

16 


(17) 


The  value  C.  is  a  distance  in  units  of  nautical  miles 
i 

or  yards  in  each  of  the  cases  above. 


hi 


The  conversion  to  ground  speed  was  easily  accomplished;  (l) 
for  C.  in  nautical  miles  and  (2)  for  C.  in  yards; 


Ci(nm)   3600(sec)    3600-Cj. 

777^7 V  ~ TTT ^ ~     =  t-  .  ,,  r~ knot 


(1)  GS  "  TIME(sec)-   l(hr) 


TIME 


(2)  GS  = 


Ci  (yds)  3600(sec)    l(nm)    =  18-C-, 
TIME (sec)    1 (hr)    2000(yds)  '  TIME 


(18) 


s  (19) 


where  TIME  is  the  time  interval  in  seconds  between  each 
range  measurement. 
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VII.   TRUE  COURSE  DETERMINATION 

The  true  course  determination  used  the  ground  speed  as 
calculated  in  the  previous  section.  Figures  6  and  7  of  Sec 
tion  V  provided  the  geometry  for  the  problem  and  Figure  7~b 
is  annotated  for  the  following  mathematical  development. 
The  coordinates  of  X  and  Y  were  first  determined.  This  was 
done  by  the  simultaneous  solution  of  a  pair  of  quadratic 
equations. 


(X-a)2  +  (Y-b)2  =  G2 
(X-c)2  +  (Y-d)2  =  G2 


Expanding  these  equations; 

X2  -  2aX  +  a2  +  Y2  -  2bY  +  b2  =  G2      (a) 
X2  -  2cX  +  C2  +  Y2  -  2dY  +  d2  =  G2      (b) 


Subtract  (b)  from  (a); 


2(c-a)X  +  a2  -  c2  +  2(d-b)Y  +  b2  -  d2 


=  G7.   -  G 


Solve    for    Y     in    terms    of    X; 

2(d-b)Y    =    2(a-c)X    +    G 
(  \  G 

Y  -  Td^bT  x 


-    G!       +    c2    +    d2    -    a2     -    b2 


-    G 


2       .        .2  2  .2 

+    c       +d       -a       -b 


yTd^bT 


Let 


D  a        ~       C  AC  N 

R    =    -, r-         and  S    =    


d    -    b 


2  2  2  2 

G         +    c       +    d       -    a 


-    b 


2(d-b) 


Then 


Y    =    R-X    +    S 


(20) 


kk 


Substitute  the  value  for  Y  in  equation  i(b)  above; 

X2  -  2cX  +  c2  +  (R-X  +  S)2  -  2d(R-X  +  S)  +  d2  =  G^ 
X2  -  2cX  +  c2  +  R2X2  +  2RSX  +  S2  -  2dRX  -  2dS 

2        2 

+  d   -  G£   =0 
Col  1 ect  i  ng  terms  ; 

2    2  2    2    2  2 

(1+R  )X   +  (2RS-2c-2dR)X  +  (c  +d  +S  -2dS-GE  )  =  0 
This  is  a  quadratic  in  X; 


X  = 


-B  +   B  -kAC 
2A 


(21) 


where 


A  =  (1+R  ) 

B  =  (2RS-2c-2dR) 

,222  2  v 

C  =  (c  +d  +S  -2dS-G£  ) 

Equation  (18)  is  a  quadratic  form  providing  two  possible 

solutions.   The  values  of  the  coordinates  X  and  Y  determined 

by  both  solutions  of  X  provide  points  which  are  mirror  images 

on  each  side  of  a  line  joining  the  points  (a,b)  and  (c,d) 

in  Figure  7~b.   The  choice  of  coordinates  (X,Y)  depends  on 

the  initial  heading  of  the  aircraft  and  whether  a  right  or 
left  turn  is  executed  on  the  second  leg. 

From  the  coordinates  (X,Y)  the  drift  vector  was  then 
determined.   The  coordinates  (e,f)  are  known. 

(22) 


W  =  ^j(X-e)2  +  (Y-f) 


The  drift  vector  was  then  applied  to  the  terminal  heading 
to  determine  the  true  course  of  the  aircraft  for_.  accurate 
track  determination. 
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VIII.   CONCLUS IONS 


The  study  showed  that  the  method  of  determining  sono- 
buoy  location  relative  to  the  aircraft  with  consecutive, 
s hor t - i n te rva 1  range  measurement  techniques  is  feasible. 
Moreover  the  application  of  the  principle  of  the  first 
derivative  of  the  range  with  respect  to  time  of  a  buoy  near 
the  aircraft  heading  as  a  means  for  determining  ground  speed 
is  accurate  and  relatively  insensitive  to  range  errors,  num- 
ber of  ranges  used  and  the  relative  angle  of  the  buoy. 

An  increase  in  the  number  of  consecutive  ranges  mea- 
sured in  excess  of  three,  reduces  the  error  in  buoy  location 
but  the  reduction  is  slight.   The  aircraft  should  not  sac- 
rifice flexibility  to  acquire  additional  range  measurements 
unless  tactics  permit. 

In  the  sensitivity  analysis  of  the  sonobuoy  location 
technique,  the  convergence  process  appeared  to  break  down 
for  acute  relative  buoy  angles  where  the  estimate  of  the 
relative  angle  lay  between  the  actual  relative  angle  and 
the  aircraft  heading.   Before  implementation  of  this  system 
further  research  should  be  conducted  in  this  area  to  pos- 
sibly provide  a  3&0   zone  of  reliable  information  on  buoy 
location.   The  model,  at  present,  can  provide  3&0   of  re- 
liable information.   For  acute  relative  buoy  angles,  with 
respect  to  the  aircraft  heading,  however,  the  estimated 
angle  must  be  greater  than  the  actual  angle. 

It  is  felt  that  the  utilization  of  buoy  ranging  in-, 
formation  for  the  location  of  sonobuoys,  apart  from 


A6 


t r i angu 1  a t i on  procedures,  has  a  broad  field  of  application, 
yet  uncovered.   A  preliminary  study  was  conducted  to  deter- 
mine relative  buoy  positions  independent  of  aircraft  navi- 
gation.  Ranges  to  all  buoys  in  the  pattern  were  measured 

from  a  single  aircraft  position.   For  a  pattern  consisting 

N! 


of  N  buoys  there  were 


possible  two-buoy  comb i - 


(21)  (N-2)  ! 
nations.   The  actual  buoy  positions  were  not  known  thus  the 

assumed  positions  were  used.   The  ranges  were  plotted  in 

reverse  from  the  assumed  positions  of  the  buoys  generating 

N  ! 
aircraft  fix  positions  determined  by  the  ranges 

(2!)  (N-2)  ! 

from  each  two-buoy  combination.   From  these  resulting  fixes 

an  average  position  was  determined  by  least  squares  method. 

An  iterative  process  was  then  developed  whereby  the  assumed 

buoy  positions   were  adjusted  such  that  the  reverse  plot- 

N! 


ting  of  the  ranges  would.generate 


aircraft  fixes 


(21)  (N-2)  1 
that  were  coincident.   This  did  not  produce  a  unique  set 

of  final  assumed  positions  for  the  buoys  but  did  establish 
a  locus  for  each  buoy  in  the  pattern  relative  to  the  air- 
craft . pos i t i on .   It  is  felt  that  further  study  in  this  area 
would  provide,  possibly,  a  more  versatile  use  of  sonobuoy 
ranging. 
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APPENDIX  A,  SONOBUOY  LOCATION  ERROR  ANALYSIS 

The  first  study  conducted  with  respect  to  sonobuoy 
location  was  to  determine  the  response  of  the  speed  of 
convergence  of  the  iterative  process  as  affected  by  the  var- 
iance of  ; 

(1)  The  relative  angle  of  the  buoy  with  respect  to  the 
aircraft  heading,  6. 

(2)  Varying  estimates,  G  ,  of  the  relative  angle. 

(3)  Range  from  the  aircraft  to  the  buoy. 

Figures  8  through  11  show  the  effect  on  the  number  of  itera- 
tions required  to  converge  by  variance  in  the  estimate  of 
the  relative  angle  of  the  buoy.   The  actual  relative  angle 
ranged  from  10   to  80   in  ten  degree  increments.   The  esti- 
mate of  the  relative  angle,  in  all  but  one  case,  ranged 
from  -15   to  +15   in  error  of  the  actual  relative  angle  in 
increments  of  one  degree.   For  the  case  where  the  relative 
buoy  angle  equals  10   the  angle  estimate,  0  ,  started  at  0 
which  was  -10   in  error  of  the  actual  angle.   The  convergence 
is  quite  rapid  for  all  estimates,  0  ,  particularly  when 
6  m  6 .   The  convergence  process  breaks  down  when  the  buoy 

is  close  to  the  heading  of  the  aircraft  and  an  estimate  8  , 

a  o 

between  the  actual  relative  angle  and  the  aircraft  heading 
is  chosen.   This  can  be  seen  in  Figure  8.   It  should  be 
noted  that  for  any  relative  angle,  8,  however  small,  if  an 
estimate,  8  ,  is  chosen  which  is  less  acute  than  q,  the 
process  provides  rapid  convergence.   Data  for  the  effect 
on  the  speed  of  convergence  for  variance  in  buoy  range  is 
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not  included.   The  effect  was  relatively  insignificant. 
Buoy  ranges  at  30  nm.  and  5  nm.  were  tested  with  at  most 
one'  iteration  difference  (in  favor  of  the  30  nm.  range)  in 
the  convergence.   The  data  in  Figures  8  through  11  repre- 
sent a  buoy  range  of  approximately  15  nm.  and  a  range  error 

of  zero.   The  convergence  process  stopped  when  the  buoy 

. . — — 

location  error  was  within  .5%  of  the  range  to  the  buoy. 

The  second  study  incorporated  a  range  error  to  deter- 
mine its  effect  on  final  buoy  location   at  the  termination 
of  the  convergence  process.   Relative  buoy  angles  of  20  , 


^5   and  80   were  investigated.   Figures  12,  13  and  \k    sho 
the  total  buoy  location  error  in  yards  versus  the  appli- 
cation of  the  range  error  in  the  following  manner  to  the 
th  ree  ranges  ; 

R    R    R 


w 


0  0  0  no  error  applied 

+  +  +  all  errors  added 

-  all  errors  subtracted 
+  -  +  other  combinations 

-  +  -  other  combinations 

There  appears  to  be  an  approximate  one-to-one  relationship 
of  the  average  buoy  location  error  to  the  range  error  ap- 
plied which  is  appealing  with  respect  to  the  small  angles 
utilized  between  range  measurements. 
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FIGURE   8 

NUMBER  OF  ITERATIONS  REQUIRED  FOR  CONVERGENCE  VERSUS 
ESTIMATE  OF  RELATIVE  BUOY  ANGLE  THETA 
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APPENDIX  B,  GROUND  SPEED  DETERMINATION  ERROR  ANALYSIS 

The  results  of  the  study  conducted  on  the  ground  speed 
determination  are  pictured  in  Figures  15  through  20.   They 
represent  the  ground  speed  error  in  knots  versus  the  range 
error  applied,  under  the  varying  conditions  listed  below; 

(1)  Number  of  ranges  measured  (3  ,  A  ,  .  .  .  ,  8)  . 

(2)  Time  between  range  measurements. 

(3)  Relative  bearing  of  buoy  to  aircraft  heading. 
In  no  case  does  the  ground  speed  exceed  1.75  knots. 

The  case  of-the  relative  buoy  angle  of  zero  was  not  tested. 
It  was  felt  that  this  was  a  trivial  case.   The  change  in 
horizontal  range  with  respect  to  time  of  a  buoy  dead  ahead 
was  the  ground  speed  exactly.   It  was  desired  to  examine 
the  accuracy  of  the  method  of  application  of  first  deriva- 
tive of  range  with  respect  to  time  as  the  relative  buoy 
angle  increased  from  zero. 
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FIGURE  16 

GROUND  SPEED  ERROR  IN  KNOTS  VERSUS  ERROR  IN  RANGE  DUE  TO 
DRIFT  IN  BUOY  FREQUENCY  WITH  RESPECT  TO  TIME  -  k    RANGES 
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FIGURE  17 

GROUND  SPEED  ERROR  IN  KNOTS  VERSUS  ERROR  IN  RANGE  DUE  TO 
DRIFT  IN  BUOY  FREQUENCY  WITH  RESPECT  TO  TIME  -   5  RANGES 
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FIGURE  18 

GROUND  SPEED  ERROR  IN  KNOTS  VERSUS  ERROR  IN  RANGE  DUE  TO 
DRIFT  IN  BUOY  FREQUENCY  WITH  RESPECT  TO  TIME  -   6  RANGES 
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FIGURE  19 

GROUND  SPEED  ERROR  IN  KNOTS  VERSUS  ERROR  IN  RANGE  DUE  TO 
DRIFT  IN  BUOY  FREQUENCY  WITH  RESPECT  TO  TIME  -   7  RANGES 
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FIGURE  20 

GROUND  SPEED  ERROR  IN  KNOTS  VERSUS  ERROR  IN  RANGE  DUE' TO 
DRIFT  IN  BUOY  FREQUENCY  WITH  RESPECT  TO  TIME  -   8  RANGES 


63 


SONOBUOY  LOCATION  PROGRAM 


START 


INPUT 

1  Estimate  of  buoy  relative 
angle  ,  0Q 

2  A/C  coordinates,  A.,B. 

3  Measured  ranges,  R 


i   i 


I 


Calculate  initial  buoy 

coord  i  na  tes  X  , Y 
o   o 

1  n  1  n 

:   =  -  y  R.cos6n  -  -  7  A 
o    n  .  L  .,  i     o    n  .  L  , 


i=l 

n 


o    n  .  l. 


i=1 

n 


i=l 


.  s  i  no 
i     o 


n  .  L « 


LOOP 
J  iterations 


Calculate  buoy  coordinates  X  &  Y. 
X  and  Y  =  polynomials  in  terms  of 
X  and  Y.   Loop  generates  iterative 
convergence  process  using  previous 
values  of  X  and  Y  to  calculate  new 
values. 


no 


-  - 

\ 

yes 

1 

OUTPUT 
Buoy  1  oca  t  i  on  s 
X,T 

' 

1 

STOP 

6^ 


SONOBUOY  LOCATION  SIMULATION  PROGRAM 


START 


INPUT 
1  A/C  coordinates,  (A . B . ) 


i  .  i 


2  Actual  buoy  coordinates,  (C,D) 

3  Estimate  of  relative  buoy 
angle,  0^ 


Calculate  actual 
ranges  to  buoy 

Rj  =  ^(C-A.)2+  (D-B.) 


I 


Establish  acceptable  error 
test  value  as  .5%    of  range 

TESTX  =  .005(R.) 

TESTY  =  .005(R.) 


i=1     i  =  1     i=1  ' 


Calculate  initial  values  of 

buoy  coordinates  using  9  Q 

1  n  n 

XQ  =  -  I    R.cosO^  -Ma. 

°    n  .  L  ,  i     o  n  .  ^  ,  i 


=  1 

n 


i  =  l 

n 


o    n .L .     i     o    n  .L , 
i=l  i=1 


LOOP 
tera  t  i  ons 


I 


Calculate  buoy  coordinates  X  &  Y. 
X  and  Y  =  polynomials  in  term-s  of 
X  and  Y.   Loop  generates  iterative 
convergence  process  using  previous 
values  of  X  and  Y  to  calculate  new 
values 


0 


no 


no 


Establish  error  in 
each  iteration 

ERRORX  =   C-X 
ERRORY  =   D-Y 


OUTPUT 
X,Y  values  of  buoy  coordinates 

and  number  of  iterations 
required  to  converge  to  within 
acceptable  values 


STOP 
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GROUND  SPEED  DETERMINATION  PROGRAM 


START 


INPUT 
Measured  ranges,  R. 


N  =  number  of  ranges 


no 


C  =' 


2      2    2 

R  -2R  +R, 

1      2    3 


2  2  2  2  2  2 

2R. -2R    -R ,+R-R    +  R 

1  2  3  it  5  6 


C    =' 


R1-R2  +  R3-Rlt-R5  +  R26-R27  +  R 


IF 


Convert  distance  C  to 
g  round  speed 

18xC 


GNDSPED  = 


TIME 


OUTPUT 

G  round  speed 

i  n  knot  s 


STOP 
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SONOBUOY  LOCATION  PROGRAMo   WI TH .MOD  I FIC ATI ONS  ,  TO  BE  USFD 
FOR  ACTUAL  DETERMINATION  OF  BUOY  COORDINATESo   INPUT  VALUES 
ARE  ACTUAL  RANGE  MEASUREMENTS  TAKEN,  THE  ASSOCIATED  A/C 
COORDINATES  AT  TIME  RANGES  TAKEN  AND  AN  ESTIMATE  OF  THE  REL- 
ATIVE BUOY  ANGLE  WITH  RESPECT  TO  A/C  HEADINGo 


l 

IMPLICIT  REAL*8  (A-H,0-Z) 

REAL*8    DSQRT,DCOS,DSIN,DATAN,PUO),PPUO),QUO)  , 
10QUO)  tSUC)  »  SSUO)  tTUC)  »R(10  )  ,  A(I0  ),  BUG  ),  DABS 

1  FORMAT(FlOoA) 

2  FORMATUOF80A) 

3  FORMATUOP80A) 

4  F0RMATU0F10o5) 

5  FORMAT  (,0,T7,«  X»,T17,«Y',T24,»  ERRORX1  ,T34,  'ERRORY'/) 

6  FORMATU     »4F10o4) 

INPUT  THE  ESTIMATE  OF  THE  RELATIVE  ANGLE  OF  THE  BUOYo 

READ(5,l)ThETA 
N  =  3 
Z  =  N 

INPUT  THE  A/C  COORDINATES  (A,B>  AT  EACH  RANGE  MEASUREMENTo 

READ(5,2)( A( I ) ,1=1, N) 
READ(5,3) (B(I),I=1,N) 
READ<5,4) (R(I) ,1=1, N) 
ABAR=Oo 
BBAR=Oo 
RBAR=Oo 
DC  50  1=1, N 
ABAR=ABAR+AU  ) 
BBAR=BBAR+B( I ) 
RBAR=RBAR+R< I ) 
50  CONTINUE 

CALCULATE  INITIAL  VALUES  OF  X  AND  Y  (BUOY  COORDINATES)  USING 
ESTIMATED  VALUE  OF  THE  RELATIVE  BUOY  ANGLEo 

RAD=57o2^577951 
THETAR=THETA/RAO 

X=(RBAF*DCOS(THETAR)+ABAR) /Z 
Y=(REAR*DS  IMTHETAR  )  +  BBAR)/Z 
WRITE(6,5) 

COMMENCE  ITERATIVE  PROCESS  FCR  CONVERGENCE  ON  VALUES  FOR  THE 
BUOY  COORDINATESo   NUMBER  OF  ITERATIONS  MAY  BE  INPUT  HEREo 

DO  100  J=l,15 
-  Vl=Oo 
V2=0o 
V3  =  0o 
V4  =  0o 
V5  =  0o 
V6  =  0o 
V7  =  0o 

DO  70  1=1, N 

P(  I  )  =  (  (X-A(  I  )  )**2+(Y-B( I)  )**2)**lo5 
PPU  )=(  (X-MII  )**2+(Y-B(  I)  )**2)**o5 
Q(I)=(X-A( I ) )**2 
OC( I )=X-A(  I  ) 
S(I)=(Y-B< I ) )**2 
SSU)  =  Y-B(I) 
T (I)=(X-A(  I )  )+(Y-B<  I  )  ) 
V1  =  V1  +  R(  I  )*S(  I  )/PU  ) 
V2=V2+R(I  )*T<  I  )/P(  I  ) 
V3  =  V3+R(  I  )*CU  )/P(  I  ) 
VA=V4+R( I )*G0( I )/PP(I ) 
V5=V5+P( I )*SS< I )/PP(I ) 
V6  =  V6+AU) 
V7  =  V7  +  BU  ) 
70  CONTINUE 
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100 


D1=1-V1/Z 

D2=V2/Z 

D3  =  D2 

D4=l-V3/Z 

D5=(Y*V2-X*VH-V4+V6)/Z 

D6=(X*V2-Y*V3+V5+V7)/Z 

X=(D5*D4-D3*D6)/(01*D4-D3*D2) 

Y=(D1*D6-D5*D2)/(D1*D4-D3*D2) 

CCNT  INUE 


WRITE 
ATIV6 


OUT    FINAL 
PROCESSo 


VALUES    OF    BUOY    COORDINATES    AT    END    OF    ITER- 


WRITE(6T6)XTY 
STOP 

END 
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SOMOBUOY  LOCATION  SENSITIVITY  ANALYSIS  PROGRAM,,   DATA  OUTPUT 
FROM  PROGRAM  DETERMINES  THF  RATE  OF  CONVERGENCE  TO  WITHIN 
ACCEPTABLE  VALUES  OF  THE  SONOBUOY  COORDINATES  AS  AFFECTED  BY 
VARIANCE  OF  THE  FOLLOWING*  NUMBER  OF  RANGE  MEASUREMENTS  USED 
ANC  THE  AVERAGE  R/5NGE  TO  THE  BUOY,  RELATIVE  ANGLE  OF  THE 
SONCBUOY  TC  THE  AIRCRAFT  HEADING  AND  RELATIVE  BEARING  ESTI- 
MATEo 

IMPLICIT  RE£L*8  (A-HtOZI 

REAL*8  DSQRT,DCOS,DSIN?DATAN,P(10) , PP( 10 ) , C ( 10 ) , 
lQQ(lO), S(10)tSS(10l ,T(10),R(10) ,A( 10) ,B< 10) 

1  F0RMAT(3F10o4) 

2  F0RMAT(10F8o4) 

3  F0RMAT(10F8o4) 

5  FORMATS  *0' ////) 

6  FORMAT ( '0»T5,  'THETA « , T17 , « X* ,127,"^ ,T37,«C« ,T47,«D»/) 

7  FORMAT! • 0« 5F10o4//) 

8  FORMAT! •     «2F10o4,I10) 

ESTABLISH  THE  NUMBER  OF  RANGE  MEASUREMENTS  TAKENo 

N  =  7 

Z  =  N 

REAC  INPUT  VALUES  OF  A/C  COORDINATES  FOR  EACH  OF  THE  RANGE 
MEASUREMENTSo 

READ(5t2)( A(I) ,1=1, N) 
READ(5,3H B(I ),I=1,N) 
GO  TO  400 
200  WRITE(6,5) 

REAC  INPUT  VALUES  OF  THE  RELATIVE  BEARING  ESTIMATE  AND  THE 
ACTUAL  BUOY  COORDINATES*) 

400  READ(5,1)THETA,C,D 

IF(CoEQolOCO)GO  TO  300 

ABAR=0 

BBAR=0 

DETERMINE    ACTUAL    RANGES    FROM    A/C    "10    BUOYo 

DO  50  1=1, N 

R(I)=DSQRT( (C-A( I ) )**2+(D-B( I)  )**2) 
ABAR=ABAR+A< I > 
BBAR=BBAR+B(I ) 
50  CONTINUE 

CALCULATE  ACCEPTABLE  ERROR  (TEST)  AS  o5?  OF  RANGE  TO  BUOYo 

TESTX=R(2)/200o 
TESTY=TESTX 
RAD=5702<;577951 
RBAR=0 
DO  60  1=1, N 
RBAR=RBAR+R(I ) 
60  CCNTINUE 

DETERMINE  INITIAL  VALUES  OF  BUOY  COORDINATES  USING  ESTIMATED 
VALUE  CF  RELATIVE  BEARINGo 

THETAR=THETA/RAD 
X=(RBAR*DCCS(THETAR )+ABAR)/Z 
Y=(RBAR*DSIMTHETAR)  +  BBAR) /Z 
WRITE(6,6) 
WRITE(6,7)THETA,X,Y,C,D 
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commence  iterative  process  to  converge  on  actual  values  of 
bucy  coordinate-So 

DO  100  J=l,20 

V1=0 

V2  =  0 

V3  =  0 

VA=0 

V5  =  0 

V6  =  0 

V7  =  0 

CALCULATE  COEFFICIENT  MATRIX  AND  CONSTANTS  FOR  SYSTEM  OF 
NORMAL  EQUATIONS  USED  IN  THE  ITERATIVE  PROCESSo 

DC    70    1=1, N 

P(II=((X-A(I)I**2+(Y-B(II )**2)**lo5 
PP(I)=( (X-A(I)  )**2+<Y-BU  )  l*:*2)^o5 
Q(I  )=<X-A<  I  )  )**2 
QQU  )=X-A(  I  ) 
S(I)=(Y-B(  I  )  )**2 
SS(I)=Y-B(I) 
T(I )=(X-A< I ) )*(Y-B< I ) ) 
V1=V1+R( I )*S( I )/P( I ) 
V2  =  V2  +  R( I)*T( I  )/P( I) 
V3=V3+R( I )*C( I )/P< I  ) 
V4=V4+R( I )*CQ( I )/PP(I ) 
V5=V5+R< I)*SS< I )/PP(I ) 
V6=V6+A< I > 
V7=V7+B(I ) 
70    CONTINUE 
Dl=    1-V1/Z 
D2=V2/Z 
D3  =  D2 
D4=l-V3/Z 

D5=< Y*V2-X*V1+V4+V6)/Z 
D6=<X*V2-Y*V3+V^+V7)/Z 

USE  KRAMER'S  RULE  WITH  RESPECT  TO  THE  COEFFICIENTS  CALCULAT- 
ED ABOVE  TO  DETERMINE  VALUES  OF  BUOY  COORDINATESo 

X=(D5*D4-D3*D6)/(D1*D4-D3*D2) 

Y=(D1*D6-D5*D2)/(D1*D4-D3*D2) 
WRITE(6,8)XtYt J 

IF  ERROR  IS  WITHIN  ACCEPTABLE  VALUE,  STOP  ITERATIVE  PROCESS 
AND  NOTE  THE  NUMBER  OF  ITERATIONS  REQUIRED  TO  CONVERGE  AND 
THE  FINAL  BUOY  COCRDINATESo 

~ TESTXX=DABS(C-X) 

TESTYY=DABS(D-Y) 

IF(TESTXXoGToTESTX)G0    TO    100 

IF(TESTYYoGToTESTY)G0    TO    100 

GO  TO  200 
100    CONTINUE 

GO  TO  200 
300    STOP 

END 
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GROUND  SPEED  DETERMINATION  PROGRAM^   FOR  USE  IN  THE  ACTUAL 
DETERMINATION  OF  GROUND  SPEED  OF  THE  AIRCRAFTo   INPUTS 
SUPPLIED  ARE  TIME  BETWEFN  RANGE  MEASUREMENTS  AND  THE 
ACTUAL  RANGES  MEASUREDo 

IMPLICIT  REAL*8  (A-H,0-Z) 
REAL+8  DSQPTtR(lO) 

1  FORMAT(10F8o4) 

2  FORMAT! ■0«F12o4) 

ESTABLISH    NUMBER    OF    RANGES    TAKEN    AND    TIME    BETWEEN    RANGES     IN 
SECCNDSo 

N  =  3 
Z  =  N 
TIME=3o 

INPUT  MEASURED  RANGESo 

READ(5»1)(R(I ) tI=l,N) 

IF(NoGTo3)GC    TO    8 

GS=DSQPT((R(1 )**2-2o*R(2)**2+R(3)**2)/2o  ) 

GO    TO    20 

8  IF(N0GTo4)GC    TO    9 

GS=DSQRT((R(1 )**2-R ( 2)**2-R< 3 )**2+R< 4 )**2 ) /43 ) 
GO    TO    20 

9  IF(NoGTo5)GC    TO    10 

GS  =  DSQRT((2o*Rm**2-3o*R(2)**2  +  R(3)**2-R(4)**2  + 
lR(5)**2)/8o  ) 
GO    TO    20 

10  IF(NoGT06)GC    TO    11 

GS=DSQRT(( 2o*R(l )**2-2o*R( 2)**2-P ( 3)**2+R < 4)**2 
1-R(5)**2+R(6)**2)/12> ) 
GO    TO    20 

11  IF(NoGT07)GC    TO    12 
GS=DS0RT<<R<l)**2-R(2)**2+R(3)**2-2o*R(4l**2+R(5)**2 

1-R(6)**2+R<7)**2  )/12->  ) 
GO   TO    20 

12  GS-DSQPT((R(1  )**2-R  (  2  )  **  2+R  (  3  )**2-R  (  4  )**2-R<  5  )**2 
1+R(6)**2-R< 7)**2+R(8)**2)/160 ) 

20    GN0SPD=GS*18o/TIME 
WRITE(6,2)GNDSPD 
STOP 
END 


73 


GROUND  SPECD  DETERMINATION  SENSITIVITY  ANALYSIS  PROGRAM^ 
CATA  OUTPUT  RELATES  TO  GROUND  SPEED  ERROR  AS  AFFECTED  RY 
VARIANCE  OF  RANGE  ERRORS,  RELATIVE  BUOY  ANGLE,  RANGE  TO 
BUOY,  NUMBER  OF  RANGES  MEASURED  AND  TIME  BETWEEN  RANGE 

MEASUREMENTSo 

i 
i 

IMPLICIT    REAf*8     (A-H,0-Z) 

REALMS    DSQRT,DABS,RUO)  ,A(  10  ) ,  B  ( 10  )  ,  ERR  (  10  ) 

1  FORMAT(2F10o4) 

2  F0RMAT(8F10o4) 

5  FORMATS 'O'Tl , • GROUNDSP EED« T 14,  'ERROR'  ) 

6  FORMAT( •  «2F10oA) 

ESTABLISH  ACTUAL  DISTANCE  FLOWN  IN  NAUTICAL  MILES  BETWEEN 
RANGE  MEASUREMENTS,  TIME  BETWEEN  RANGE  MEASUREMENTS  AND 
ACTUAL  GROUND  SPEED  TO  BE  USED  AS  A  TEST  FOR  DATA  GROUND- 
SPEEDo 

E  =  o25 

TIME=4o5 

TEST=E*3600o/TIME 

N  =  6 

Z  =  N 

INPUT    ACTUAL    BUOY    COORDINATES    TO    SIMULATE    VARYING    BUOY 
ANGLES    AND    RANGES    RELATIVE    TO    THE    AIRCRAFTo 

READ(5,1 )C,D 

INPUT  AIRCRAFT  COORDINATES  AND  RANGE  ERROR  APPLIED  TO 
ACTUAL  RANGE  TO  SIMULATE  MEASURED  RANGE  WITH  ERPORo 

READ(5,2)(  Ad  )  ,1=1, N) 

READ(5,2H  B(I  )  ,1=1, N) 

WRITE(6,5) 

DO  100  J =1,70 

READ(5,2)(ERR( I) ,1=1, N) 

CALCULATE    MEASURED    RANGE"  AS    DIFFERING    FROM    ACTUAL    RANGE 
BY    ERROR    APPLIEDo 

DO    200    1=1  ,N 

R(I  )=DSQRT(  (OA(  I)  )**2  +  <D-B(  I)  )**2) 

RU  )=R(  I  KERR(I) 
200    CONTINUE 

CALCULATE    VALUE    OF    GROUND    SPEED    FROM    INPUT    DATAo 

IF(NoGTo?)GC    TO    30 
_  GS=DSQRT(( R(l»**2-2o*R(2)**2+R<3)**2)/2o ) 
GO    TO    50 

30  IF(NoGTo4)G0    TO    31 

GS=DSQRT( (R(l)**2-R(2)**2-R(3)**2+R(4)**2)/4o ) 
GO    TO    50 

31  IF(NoGT05)GC    TO    32 

GS=DSOPT( (2o*R(l)**2-3o*R(2)**2+R( 3)**2-R ( 4)**2 
1  +  R(5  )**2)/8o ) 
GO    TO    50 

32  IF(NoGTo6)GC    TO    33 

GS=DSORT( ( 2o*R(l)**2-2o*R(2)**2-R(3)**2+R<4)**2 
1-R(5)**2+R(6)**2)/12, ) 
GO   TO    50 

33  IF(NoGTo7)G0    TO    34 

GS  =  DSQRT( ( R(l  )**2-R ( 2) **2+R( 3)**2-2- *R <4)**2+R< 5 )**2 
l-R(6)**2+R(7)**2)/12o  > 
GO    TO    50 

34  GS  =  D SOR T  (  (  R  (  1  )  **  2- R  (  2  )  **  2+  R  (  3 )  **  2- R  (  4  )  **  2- P  (  5  )  **2 
l+R(6)**2-R(7)**2+R(8)**2)/16o ) 

50    GNDSPD=GS*3600o/TIME 


Ik 


TFST  CALCULATED  GROUND  SPEED  AGAINST  ACTUAL  GROUND  SPEED  TO 
DETERMINE  ERRORo 

ERROR=DABS( TEST-GNDSPD) 
100  WRITE(6,6)GNDSPD, ERROR 
STOP 
END 
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